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The search for cosmic strings has been of renewed interest with the advent of precision cosmology. 
In this note we give a quantitative description of the nonlinear matter density fluctuations that can 
form from a scaling network of cosmic string wakes. Specifically, we compute the distribution of 
dark matter halos. These halos would possess strong correlations in position space that should have 
survived until today. We also discuss the challenges involved in their detection due to their small 
size and the complex dynamics of their formation. 
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I. INTRODUCTION 

Cosmic strings are linear topological defects predicted 
by many models of particle physics that go beyond the 
Standard Model, e.g. in supergravity [1:, brane inflation 
[5] and "String Gas Cosmology" [3] models. In models 
which admit these defects, a network of cosmic strings 
is inevitably [4] formed during a phase transition in the 
very early universe and will persist to the present time. 
This is true as long as the phase transition occurred af- 
ter the period of inflation (provided there was inflation) . 
Since strings carry energy, they lead to gravitational ef- 
fects which induce distinctive signatures in cosmological 
observations. Since the string tension (which equals the 
mass per unit length fi) increases as the energy scale r] 
at the time of string formation grows, the magnitude of 
the predicted signatures of strings increases as the energy 
scale grows. Thus, a search for observational signatures 
of strings allows us to constrain the physics occurring at 
energies much larger that can be probed by earth-based 
particle accelerators. 

It has long been known that cosmic strings are able to 
produce many interesting effects for cosmology. In par- 
ticular, cosmic strings lead to a scale-invariant spectrum 
of cosmological perturbations 5 . Initially, the focus of 
interest was on strings which have large enough tension 
to explain the entire amplitude of the power spectrum 
[5] . Such strings would need to have a value of the ten- 
sion given (in dimensionless units) by G/i ~ 10~ 6 , where 
G is Newton's gravitational constant. The perturbations 
produced by cosmic strings are, however, active and in- 
coherent [5] and hence do not lead to acoustic oscilla- 
tions in the angular power spectrum of cosmic microwave 
(CMB) anisotropies. Thus, when conclusive evidence for 
the existence of these oscillations was reached [7J, it be- 
came clear that cosmic strings could not be the dominant 
source of structure formation 8 and interest in cosmolog- 
ical effects of strings collapsed. However, in light of the 
fact that in many models of inflation (or other explana- 



tions for the dominant Gaussian component of the power 
spectrum of cosmological perturbations) cosmic strings 
are predicted, there has lately been a revival of interest 
in searching for cosmic strings (see e.g. [5] for a recent 
survey). A second reason for the revival of interest in 
searching for strings is that there have been great ad- 
vances and technological breakthroughs in observational 
cosmology which are now making it possible to search for 
signals of cosmic strings with much lower tensions than 
those previously considered. 

The searches have been progressing in many directions. 
First of all, cosmic strings will contribute to the power 
spectrum of density fluctuations [5] and to associated 
CMB anisotropies [TU]. Signatures of these fluctuations 
can be searched for in optical and infrared galaxy sur- 
veys and in CMB temperature anisotropy maps. Cos- 
mic strings leave behind distinctive signals in 21cm red- 
shift surveys [TT] and in CMB polarization maps |12j . 
Considerable work has been done to study the effects of 
strings in the angular power spectrum of CMB tempera- 
ture and polarization maps [13] . Analyses combining the 
SPT [14] and WMAP 7-year [15] data were able to place 
the bound G/i < 1.7 x 10 -7 16J 1 , and a similar bound 
can be obtained [IB] using data from the ACT telescope 
[19) . Improved bounds might be achievable by analyz- 
ing CMB maps using statistical tools which are designed 
to pick out the string-induced non-Gaussianities (see e.g. 
[2TJI |2"T] for some recent studies) . 

In addition to these purely gravitational effects, cosmic 
strings can be responsible for the production of highly 
energetic bursts of particles [22], electromagnetic radia- 
tion in a wide range of frequencies [23] and they can help 
seed coherent magnetic fields on galactic scales [24 . Cos- 
mic string loops decay by emitting gravitational radiation 
[25] , Hence, strings lead to a scale-invariant stochas- 
tic background of gravitational waves [26] which can be 
constrained by pulsar timing measurements. Cusps on 
string loops may lead to additional emission of gravita- 
tional waves |27| , and a resulting constraint on the cosmic 
string tension of G/i < 4 x 10 -9 [28] has been reported. 
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1 See 1171 for earlier bounds using only the WMAP data. 
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However the bound is sensitive to details of cosmic string 
cusps which are subject to potentially large back-reaction 
effects [29j . A bound of G[i < 5.2 x 10~ 7 resulting from 
constraints on the amplitude of gravitational radiation 
from pulsar timing constraints is more realistic [30] . 

Observational optical astronomy is also experiencing 
a phase of rapid progress. Larger telescopes are prob- 
ing the universe at increasing depth, i.e. at increasing 
redshift. Since cosmic strings produce non- linear density 
perturbations at arbitrarily large redshifts, it is expected 
that the signals of cosmic strings will stand out from 
the structures produced by the main source of fluctua- 
tions (Gaussian fluctuations) more clearly at high red- 
shifts than at low ones. Signals from structures seeded 
by string loops has been studied in a recent paper by 
Schlaer et al. [31] (see also [32] for earlier work) who 
found that loops could cause significant star formation 
to occur at high redshift. Therefore cosmic strings with 
enough mass will have an impact on the epoch of reion- 
ization [33j [34] . String loops would also seed dense dark 
matter clumps whose population today could produce a 
background of gamma radiation sensitive to the Fermi 
telescope [35] . 

In this note we address early structure formation from 
cosmic string wakes. Wakes lead to cosmic structures 
with a distinctive shape in position space maps. These 
non-Gaussian features may allow the signals of strings 
to be detected in a background of Gaussian noise even if 
the amplitude of the string-induced signal is quite low, as 
studied for CMB anisotropy maps in [5D]. In Section [TT] 
we give some background about the formation and prop- 
erties of cosmic string wakes. In Section III we investigate 



the contribution of string wakes to structure formation as 
a function of redshift. Finally we discuss the implications 
and observational signals in Section [TV] 



II. BACKGROUND 

According to the Kibble mechanism [4], the formation 
of cosmic strings from a symmetry breaking phase tran- 
sition implies the existence of a string network. This 
network is expected to evolve towards a scaling solution 
which states that its properties become constant if all 
lengths are scaled to the Hubble radius [35] . The net- 
work is composed of two parts. The first is a network 
of "infinite" strings (which includes loops with radius 
greater than the Hubble radius) whose curvature radius 
R c can be shown to be of the order of the Hubble scale, 
i.e. R c — -ft with 7 ~ O(l). The fact that the curva- 
ture radius scales with t is maintained by intersections 
and self-intersections of the long strings. Such intersec- 
tions create string loops that detach from the long string 
network, thus allowing the latter to straighten out. This 
brings us to the second part of the network: a distribution 
of string loops with curvature radius much smaller than 
the Hubble scale. Following [5] we model a long string 
as a superposition of straight string segments, each of 



whose length is given by R c . This is a reasonable model 
since any wiggles on the strings would either be redshifted 
away or would cause self-intersections. The scaling so- 
lution suggests N ~ 1 long string per Hubble volume, 
which is supported by numerical simulations which yield 
N s» 1 — 10 with a loop distribution possessing a scal- 
ing peak at the loop radius of size sa 1/20 of the Hubble 
radius }37| (for earlier work on cosmic string simulations 
see e.g. [35]). 

The geometry of space around a long straight string 
will be conical with deficit angle a = 8ttG/j, [39]. A long 
string moving with speed v s will sweep a plane on which 
residing static observers would see that matter acquires a 
velocity kick v = 4:irG/iv s -/ s towards the plane (where 7 S 
is the relativistic gamma factor associated with the veloc- 
ity v s ). The streams of matter will overlap and create a 
wedge with twice the background density [30] • The over- 
dense region will eventually collapse under gravitational 
instability and form a virializcd planar structure called 
the string wake. The growth of the wake will proceed by 
matter accretion. We will briefly review the dynamics of 
this process. 

Let us consider 2 the effect of the wake on a particle's 
motion by labelling its physical distance to the wake by 
r = a(t)(x+ip(t)), where x is the comoving position, ip(t) 
the displacement and a(t) the scale factor. The initial 
conditions are given by ip(ti) — and ip(U) — —sgn(x)v. 
The Zel'dovich approximation [42] gives as equation of 
motion, 

* + ^-^ = 0. (1) 
This system can be solved to obtain, 

12 ft \ 2 / 3 12 t 2 

ip{t) = - — Gnitv B ') a ti\^—j +yG^7ru s 7 s y. (2) 

After dropping the decaying mode, we can find the dis- 
tance at which our particle decouples from the Hubble 
flow through the condition r = or more specifically 
x + 2ip = which yields, 



24 / t \ 

X ta (t) = i — GfilTVs'yati^ — J 



t \ 2/3 



(3) 



Hence w(t) = ^a(t)x ta (t), which denotes the distance 
from the wake's center of the mass shell that turns- 
around at t, is given by 



w(t) = ^-nG^v s -f s u(^^ 



t \ 4 /3 



(4) 



Matter that turned around at a height w(tt a ) will viri- 
alize at a radius \w{t ta ) showing that w(t) labels the 



Such calculations were first performed for string wakes in |41| . 
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width of the wake which contains four times the back- 
ground density of matter. Shells of infalling baryonic 
matter will collide with one another at the virial height 
and create shocks on either side of the wake as shown 
through hydrodynamical simulations [H]. The en- 
ergy of the falling particles will be thermalized in the 
wake giving a temperature 



-k B T 



1 



-mv 



shell ' 



(5) 



where v s h e ii is the speed of hydrogen particles of mass m 
at the shock. A shell that turns around at time t ta will 
hit the wake at t = (1 + 1/V2) 3 ^ 2 t ta , and therefore the 
speed will be given by 



v sh eu = r(t,t ta = t(l + l/V2)- 3/2 ) 
which yields a wake temperature of 



(G) 



T = 



16 

75 



(4(3-2V2)y^(Grf(v sls r(l) V3 (7) 



[WK](G») 2 6 (v sls 



z + l 



(8) 



with (G/x) 6 being the value of Gp in units of 10 



III. THE MASS FUNCTION 

To compute the amount of matter accreted onto wakes, 
and eventually the mass function, we must understand 
how to model the network of wakes created from our 
scaling solution. As previously mentioned, we expect N 
strings per Hubble volume. The curvature radius of the 
long strings network grows in length relative to fixed co- 
moving coordinates. This growth is realized because long 
string segments intersect and produce string loops. The 
string intercommutations occur roughly once per Hubble 
expansion time per string per Hubble volume. They also 
lead to a change of the direction of motion of a string af- 
ter intersecting another string and exchanging ends. We 
will use a toy model introduced in fS] to characterize the 
long string network analytically: we divide the time axis 
into Hubble time steps. In each time step and in each 
Hubble volume we lay down N straight string segments 
of length 7^ where 7 is less but close to unity, and each 
string has a velocity v s in a random direction. We take 
the string segments to be uncorrelated. Each string seg- 
ment lives for a Hubble time. The string network at 
different Hubble time steps is taken to be uncorrelated. 

Wakes first form behind strings at the time of equal 
matter-radiation energy density t eq , when matter pertur- 
bations can start to grow. We consider Hubble time steps 
from t eq to t labelling them t\ — t eq , t 2 — 2t eq , ...,t m = 
2 m ~ 1 t eq . At the beginning of each time step t m we create 
N string segments, at the end of the time step we annihi- 
late them. Each of these strings will create a wake with 



initial dimensions jt m x v s ^ s t m x w{t m ) (where 7 S is the 
relativistic gamma factor associated with the string ve- 
locity v s ) which will then grow in planar dimensions due 
to the Hubble expansion and in width via the matter 
accretion process discussed in the previous section. 

Consider a redshift z, then the comoving number den- 
sity of wakes laid down at time z m > z is given by 



N 



1 



#m (z m + l) 3 



0) 



Therefore the mass density in all wakes at a redshift z 
is given by summing the contribution from all the wakes 
created at different Hubble time steps: 

M w {z) = ip b (z) ^ n wake (z m )(l + z) 3 Vo\ wake (z m , z) , 

Z m >Z 

(10) 

where the last factor is the volume at redshift z of the 
nonlinear region about a wake created at a redshift z m > 
z. Since the planar dimensions of the wake expand with 
the Hubble flow, we obtain 



M w (z) = Ap b (z) ^ n wake (z m )(l + z) 3 

Z m >Z 

x {jt m x v s ^ s t m x w(z))(^ 
A — ' 45 1 + 2: 



1\ 2 



z + l 



(11) 



In Figure [I we plot the quantity M w (z)/pb(z) (where 
pb{z) is the background energy density) for the following 
parameter choices: 

Gfi = 1.5 x 10~ 7 , N = 10, v s = 1/2 and 7 = 1, (12) 

where the value of Gp is close to the current upper 
bound and the other values are taken to agree with 
those in recent cosmic string evolution simulations. This 
quantity represents the fraction of the total mass in the 
universe accreted onto wakes. 

Figure [T] shows that a considerable fraction of matter 
is accreted by wakes. Let us analyze this a bit more care- 
fully. Small density perturbations on the surface of the 
wake will grow due to gravitational instability and might 
eventually cause the wake to fragment. The dynamics 
of this process will be quite complex due to the inner 
structure of the wake: the dark matter density peaks 
at the edges while baryonic matter clusters in the cen- 
ter [43j . An analysis of the fragmentation was done in 
the simpler case of a static isothermal sheet-like cloud in 
[45] . The timescale of the perturbation growth is of the 
order of the freefall time with the length of the most un- 
stable mode equalling the sheet's width and the longest 
wavelength unstable mode being twice that. Taking the 
very optimistic route and assuming our wake behaves in 
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FIG. 1: Fraction of the total matter in the universe that is in 
cosmic string wakes at various redshifts. 



a similar way 3 , long filaments of virial diameter sa w(t) 
should form in a similar fashion to the Zel'dovich pan- 
cakes. These filaments would subsequently break into 
beads themselves, again the fastest growing mode having 
a length ss 2ww(t) [31]. Therefore the wakes will eventu- 
ally fragment into virialized halos of size ~ ir(w/2) 3 . 

With this picture in mind we can compute the comov- 
ing number density of halos of mass greater than M that 
fragments from the wakes created at redshift z m : 

n(z,z m , > M) = n wake (z m ) ^f~™ ake 0(M ha i o - M) 



){jt m x v s j s t m x w(z)) 



Zm. ~\~ 1 \ ^ 



z + l 



X 4 7 , 3 e(2 Pb (z)^w(z) 3 - M), (13) 

2 7TW yZ J o 

where O is the Heavyside step function. 

Taking into account the all wakes that exist at redshift 
z we simply get 



n(z, > M) = ^ n ( z > z ™, > M ) ■ 



(14) 



The collection of step functions given by the sum is a 
bit odd as a distribution since in reality this should 
be continuous. The discrepancy arises because of the 
model that considers the wakes to be instantaneously 
created while this process actually spans a Hubble time. 
The part of a wake that is created close to time t m 
will have different properties than the part created at 



3 The neglect of Hubble flow along the wake might not be justified. 
In an expanding universe structures with planar, cylindrical and 
spherical symmetry undergo a self-similar growth proportional 
to t A / 3 , t and t 8 / 9 respectively [46]. In this case one could argue 
that if perturbations grew as cylindrical or spherical objects, they 
would grow slower than the planar wake and therefore no clear 
fragmentation would occur. 



around 2t m which would be more similar to the begin- 
ning of the next set of wakes. Hence we should smear 
out the distribution n(z, > M) to have something real- 
istic. To do this we fix z and split the mass range in 
intervals AM m — [M(z, z m ), M(z, z m —i)] with M(z, zj) 
being the mass of the halo's formed at z from a wake 
created at the Hubble time step Zj. Since the values of 
n(z, > M{zj)) = J2 Z > z - n ( z ? z m, > M) are known, we 
can just interpolate to create a continuous n(z, > M). 
One can then obtain the comoving mass function from 

dn(z.>M) 
d log M ' 

In Figures[2]and[3j we plot and compare n(z, > M) and 
^diogM f rom our model to the standard predictions as 
found in Reed et al. [57J based on Gaussian fluctuations 
in a standard ACDM model with cosmological parame- 
ters chosen to agree with those obtained from the WMAP 
data. We use the string parameters given in Equation [T2| 
and note that the small wiggles in the plots of n(z, > M) 
are just an artefact of our interpolation. We see from the 
figures that structures formed from primordial fluctua- 
tions dominate at redshifts smaller than z ~ 20, whereas 
the string wakes dominate at higher redshifts. One thing 
to note is that the mass function for Gaussian pertur- 
bations is obtained using the extended Press-Schechter 
formalism [IS] which counts the small halos and their 
parent halos independently. For wakes, we assume here 
that smaller halos get destroyed when they merge in the 
wake to create larger ones, and hence the tail end of the 
halo distribution could be larger than what we predict 
here. 

Note that even if at some redshift the wake halos make 
up a decent fraction of the total halo number in the small 
mass region, objects of these sizes are hard to detect. 
Therefore it will be challenging to probe this mass range 
directly with any decent precision. 

The statistics studied here do not take into account 
the distinct position space correlations that wake halos 
would inherit between each other. The wake-induced ha- 
los would be concentrated in sheets, and hence might be 
visible above a Gaussian noise with significantly larger 
amplitude using topological statistics such as Minkowski 
functionals [15] ■ In the early days of interest in cosmic 
strings, the application of Minkowski functionals to dis- 
tinguish structures seeded by strings from that seeded 
by Gaussian noise was analyzed in [50] . Minkowski func- 
tionals have also recently been applied to the analysis of 
the signals of cosmic strings in 21cm redshift maps |51j . 

Note, however, that the position space correlations 
might be disrupted by the interactions between the wakes 
and the primordial perturbations, something that we 
have so far neglected. The structures corresponding 
to density fluctuations on scales larger than the largest 
wake's width cannot be affected by these interactions. 
However, the formation of these large objects might de- 
stroy any wake in their proximity. We give a rough es- 
timate of the number of surviving wakes by determining 
the fraction of space that are voids with no large halos. 
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n(>M) of DM halos at redshifts z=\ 1, 2, 3, 4) 
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n(>M) of DM halos at redshifts z={ 10, 15, 20, 25, 30} 
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FIG. 2: niz, > M) of dark matter halos for a string network 
with parameters 12 (solid) and the Reed et al. [47| predic- 
tions (dashed) at different redshifts. The colors with largest 
wavelength corresponds to the largest redshifts. 



The fraction of space that collapsed in Gaussian noise- 
induced halos of mass greater than M is given by j52j. 



F(> M, z) = 2 

Js crtt 



dS- 



1 



Erfc 



( Z) v2trt(M) 
riiiZ \) 



exp 



f"— ) 
\ 2ct(M) 2 / 



V2a(M)J' 



(15) 



Here, S cr n = 1.686(1 + z) and a(M) is the rms 
density fluctuation on a comoving length scale R = 
(3M/47rpb(0)) 1//3 which depends on the power spectrum 
of matter. Figure [4] shows F(> M w , z) for three different 
values of G/i. Here M w (z,G/j.) is the maximum mass of 
wake halos. We also plot F(> 10 3 M Q /h,z) which was 
the smallest mass that was resolved in the simulation by 
Reed et al. [47]. If F is close to 1, then most of the 
matter was accreted onto large Gaussian noise-induced 
halos, and therefore large halos wiped out any geomet- 
rical structures and spatial correlations that wake halos 
could possess. As the figure shows, the larger the red- 
shift is, the less will be the washout of string-induced 
structures by the Gaussian noise. 



Comoving Mass function of DM halo's at redshifts z={ 1, 2, 3, 4} 
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Comoving Mass function of DM halo's at redshifts z=| 10, 15, 20,25, 30} 
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FIG. 3: Mass function of dark matter halos for a string net- 
work with parameters 12 (solid) and the Reed et al.[47] predic- 
tions (dashed) at different redshifts. The colors with largest 
wavelength corresponds to the largest redshifts. 
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Fraction of matter in halos (with mass>M) at different redshifts 
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FIG. 4: Fraction of matter in halos of mass greater than wakes 
can create. We plot the curves for multiple values of G/i and 
one of a mass of 1O 3 M0 for reference. These were computed 
with the program supplied in [47] . 
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IV. DISCUSSION 

In this paper we have studied the formation of halos 
of dark matter by cosmic string wakes as a function of 
redshift. We have found that the relative importance of 
string wakes increases as the redshift increases. However, 
there are challenges in trying to detect these signals, and 
this is the topic we discuss in this section without giving 
too many details. 

The temperature of our halos can be estimated from 
the virial theorem and conservation of energy by consid- 
ering a collapsing gas cloud of radius w containing two 
times the background density pb, 



^ 2 ^virial 

3 GM 



5 r. 



virial 



3 GM 

5 ^virial 

1 2 3 GM 

7, V th ~ r 

2 5 w 



(16) 
(17) 



Here, w is the radius of our collapsing cloud [4j M — 
2pb^(w) 3 is its mass and Vth — v s h e ii is t ne initial ther- 
mal velocity in Eq. [6] The shock heated baryonic gas 
will have a temperature of 



T 



\xm v 2 

-v , 

2 k b ' 



(18) 



where m# is the mass of the proton and fi is the mean 
molecular weight [52] . Star formation can occur if the 
temperatures are high enough so that colliding molecules 
can excite each other to higher energy states and radi- 
ate away their kinetic energy. Atomic cooling requires 
T > 10 4 K while H2 cooling has channels that can be 
efficient down to temperatures as low as 200K. However 
these H2 molecules are very fragile and get destroyed by 
the radiation of the first stars, and therefore this cool- 
ing mechanism alone will not allow for multiple stars to 
form in close proximity as would be the case for halos 
formed from wakes. Figure [5] shows the temperature of 
the hottest temperatures that can be achieved with the 
parameters in Eq. 12 and /i ~ 1.22 in the case of neu- 



tral hydrogen. As we can see, the temperatures are far 
too low to have any significant star formation, and hence 
no noticeable effect on the history of reionization will be 
present. 

Gravitational lensing could also be an interesting win- 
dow considering there should be about ~ 10 6 wakes in 
our Hubble volume. However, even if each wake contains 
a lot of mass, it is spread over a large area and hence 
we would only obtain strong lensing if the background 
source lay, with us, very close to the plane delimited by a 
wake, in order to have enough mass along the integrated 
line of sight. This makes the signal too rare to be use- 
ful to statistically constrain string parameters. However, 
if one were lucky, assuming the images do not overlap, 
which could happen for extended source and a small Gfi 
or if we are not in the proximity of the wake's plane, a 



lensed source would appear as two images with no mag- 
nification or distortion. In such case one could hope to 
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FIG. 5: Redshift dependance of the temperature in the 
hottest halos that are formed from wakes. We use the pa- 
rameters in Eq. |12| and fi w 1.22. Note that in the entire 
redshift range shown, the halo temperature is higher than 
both the CMB background temperature and the average ki- 
netic temperature of the baryonic gas. 



find other similarily lensed objects along some line in the 
sky. 

Observables sensitive to the structures formed at high 
redshift would be ideal to differentiate and detect wake 
structures from the background. One such example, that 
relies on the compact objects, is gamma ray constraints 
from dense dark matter clumps [351 153 ) . However, the 
formation of these objects must proceed by a clean radial 
infall of matter in order to have a reliable density profile 
to which the gamma ray flux is sensitive. These are not 
favorable conditions in the wake since we expect strong 
tidal forces as growing perturbations oscillate around the 
center. 

In conclusion, we have shown that string wakes dom- 
inate the nonlinear structures in the universe at suffi- 
ciently high redshifts. Even at redshifts where they are 
not dominant, they may be identifiable because of the 
distinct spatial correlations which they induce. However, 
the string-induced halos are typically too small and cold 
to induce star formation, and hence it is a challenge to 
be able to find them observationally. One promising win- 
dow - studied in other works [TH ISTl I5"1HS*6"] - is via 21cm 
redshift surveys. 
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